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In the last several years, the Casimir energy for a variety of 1+1-dimensional integrable 
models has been determined from the exact S-matrix. It is shown here how to modify the 
boundary conditions to project out the lowest-energy state, which enables one to find 
excited-state energies. This is done by calculating thermodynamic expectation values of 
operators which generate discrete symmetries. This is demonstrated with a number of 
perturbed conformal field theories, including the Ising model, the three-state Potts model, 
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1. Introduction 

A wealth of physical information about a relativistic quantum particle theory is con- 
tained in the knowledge of the asymptotic particle states and their scattering matrix. This 
includes all thermodynamic quantities of an infinite- volume gas of these particles, as long 
as the gas is dilute enough so that a particle description is valid. There is a closed-form 
expression for the virial coefficients (terms in the expansion of the free energy in the fu- 
gacity) depending only on the S-matrix[|l]. Doing the calculation in this manner is usually 
unwieldy and unnecessary, because one can generally calculate thermodynamic quantities 
directly from perturbative field theory. However, the situation is markedly different in 
1+1-dimensional integrable theories. Here the S-matrix has a number of properties which 
simplify matters a great deal: it is completely elastic (momenta are conserved individually 
in a collision), the n-body S-matrix factorizes into a product of two-body ones, and one 
often can derive or conjecture it exactly 0. It is in principle possible to calculate these 
quantities using perturbed conformal field theory, but using the S-matrix is often more 
tractable, and if the S-matrix is exact, the results are non-perturbative. 

The techniques for calculating the Casimir energy (the lowest eigenvalue of the one- 
dimensional quantum Hamiltonian) were first used in the non-relativistic case in ||. Re- 
cently, there has been a great deal of interest in this field, because the study of perturbed 
conformal field theory has enabled conjectures for a large number of exact S-matrices[f|].0 
The calculation of the Casimir energy in such a case goes by the name of the Thermody- 
namic Bethe Ansatz (TBA for short) and was first done for the Yang-Lee model and the 
three-state Potts model in ref . |§ . This paper will discuss modifications of these techniques 
recently proposed in 0-0, which enable the calculation of energies of excited states. 

The result of the TBA is a set of coupled integral equations, one for each particle 
in the spectrum. These equations give the exact Casimir energy as a function of the 
volume of space. Although one can solve these equations only numerically, it is possible 
to extract either the infrared or ultraviolet asymptotics. The ultraviolet limit is where 
the theory becomes a conformal field theory, and thus the results can be compared with 
results at the conformal point. The results of this paper will be compared to conformal 
field theories with a variety of boundary conditions [fTPf 1 1 1 ] . The TBA equations have been 



studied in a large variety of exactly solvable theories, including the XXZ spin chain in a 
magnetic field[E2| and its generalizations PB| , the supersymmetric sine-Gordon model|n 



1 For a review of much of this activity, see 0. 
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all the simply-laced affine Toda theories [|15| |jT6| , a number of non-unitary mo dels P| JT7| , 



the perturbed minimal models|T^]-[j21[, and the perturbed N=l and N=2 supersymmetric 



minimal models[]22[]. Some of these results may even be experimentally verifiable []23[ . 

The TBA allows the calculation of thermodynamic quantities for a one-dimensional 
system, where the spatial dimension is a circle of circumference L, with L — > oo. Any 
thermodynamic expectation value of a <i-dimensional quantum system at temperature T 
can be alternatively viewed as an expectation value in a (d + 1) -dimensional Euclidean 
quantum field theory, where the additional dimension is a circle of circumference R = 
1/T. This is therefore equivalent to a two-dimensional theory with spacetime a torus of 
periodicity (R,L). We can view either the .R-direction or the L-direction as space, with 
the other direction then taking the role of Euclidean time. Thus the partition function 
Z(R, L) is equal to both 

Z(R,L) =tr [e~ RHL ] , (1.1) 

and 

Z(R, L) = tr [e~ LHR ] , (1.2) 

where He is the Hamiltonian for the system with space a circle of circumference C. This 
equivalence has proven very useful in conformal field theory defined on a torus, where 
it goes by the name of modular invariance, and puts severe constraints on the operators 



allowed in a theory |p4"| and on their fusion rules J25[. In the L — > oo limit, the only state 



which contributes to the partition function is the ground state. Thus 

Z{R, L -> oo) — > e~ E( - R)L , (1.3) 

where E(R) is the ground-state energy (lowest eigenvalue of Hr) of the theory on a circle 
of circumference R, and does not depend on L. In any thermodynamic system the free 
energy is F = — In Zj R, so ( |Q| ) yields 

E(R) = ^F(#,L -> oo). (1.4) 

LJ 

Therefore, the Casimir energy can be determined from the TBA. 

An important check on these results can obtained by studying the limit where the 
particle masses go to zero, so the system approaches a conformal field theory. A general 
result of conformal field theory |26|] predicts 



E{R) = ^(h + h-±), (1.5) 



where h and h are respectively the left and right conformal dimensions of the operator 
which creates the lowest energy state when acting on the vacuum. For a unitary theory 
with periodic boundary conditions (anti-periodic for any fermions present), this operator 
is the identity, so that h = h = 0. This relation gives a very useful way of checking that the 
conjectured S-matrices used in the TBA calculation are consistent — (|1.5|) relates a TBA 
result in the m — > limit to results obtainable from conformal field theory. This check 



has been successfully done in many cases||15||-[S2|, thus giving even stronger support for 



the conjectured S-matrices of these models. In addition, many of these papers numerically 
calculate corrections to (fl~5l), and compare them with results obtained with perturbation 
theory around the conformal point. 

All of this work determined E(R) for the lowest-energy state of the theory. However, 
the TBA formalism is not limited to this type of calculation. In several recent papers the 
TBA equations for a number of models were modified in order to describe the behavior 
of excited-state energies -||. To check these TBA equations, the values of h in (]1.5Q 
were calculated, and were found to match those of particular operators in the conformal 
theory. The equations were solved numerically in some cases, and the results agreed with 
the appropriate perturbed conformal field theory. 

In this paper we demonstrate two methods of deriving excited-state energies. These 
methods are general and allow the determination of excited-state energies in any model 
with a discrete symmetry. This relates the results of [§ and |J for massive theories to 



conformal results |10| , and gives a simple physical interpretation of their modified TBA 
equations. 

Both methods involve the calculation of the thermodynamic expectation value of the 
operators which implement this discrete symmetry, and are based on the way partition 
functions are calculated in conformal field theories with boundary conditions |10|] [llj. The 
unnormalized expectation value of an operator A is 

< A >= tr [Ae~ HLR ] . (1.6) 

Thus Z(R,L) =< 1 >. Just as < 1 > can be rewritten as (|TT^) by interchanging space 
and Euclidean time, we rewrite 



< A >= tr 



e -H A L 



(1.7) 



where H a is the Hamiltonian for the same model, but with new boundary conditions 
depending on the operator A. If A = 1, then the boundary conditions are periodic for 
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bosons and antiperiodic for fermions, as is standard in thermodynamics. Another simple 
example is the case of a free scalar field (p. When A is the operator which sends <fi ~^ 
Ha is the Hamiltonian for the model with antiperiodic boundary conditions. 
As in the periodic case, in the L — > oo limit 



long as the energy-momentum tensor is invariant under the action of A. However, h is not 
necessarily zero, even in a unitary theory. This is because the boundary conditions can 
project the identity state out of the spectrum. In conformal field theory, this has been 
demonstrated explicitly for the Ising model and the three-state Potts model with a variety 



There are two ways of calculating < A > when L — > oo, thus determining Ea(R)- 
The first method utilizes an imaginary chemical potential and is used when A is diagonal 
on the space of particles. The second method works when A interchanges particles, and 
involves truncating the Hilbert space. These methods can be combined for a general A. 

Section 2 is a warm-up for the rest of the paper — both methods of calculating excited- 
state energies are applied to the Ising model, where the TBA is not needed. In section 
3, the TBA in the presence of a general chemical potential is reviewed, and it is shown 
how this can be used to calculate expectation values of diagonal operators. Section 4 
demonstrates within the TBA how one truncates the Hilbert space to obtain expectation 
values of non-diagonal operators. Section 5 applies these methods to models described by 
simply-laced Toda minimal S-matrices, including the three-state Potts model. Section 6 
discusses two more models, the massless perturbation of the tricritical Ising model and a 
model which displays many of the characteristics of the minimal models. Section 7 contains 
conclusions and outlines possible future work. 

2. The Ising Model 

The field theory describing the T ^ T c Ising model in two spacetime dimensions 
consists of a single free Majorana fermion of mass m.@ We consider an ensemble of N 

2 Another field theory describing the Ising model consists of a single massive boson with S- 
matrix S = — 1. Which of the two is applicable depends on the boundary conditions, since they 
have different finite- volume thermodynamics However, the ground-state energy is the same 
in both theories, so it does not matter which theory we use here. 




of boundary conditions [|10 . 
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particles, and the rapdity 0j of the ith particle is defined by 



Pi = m sinh $i 

(2.1) 

Ei = m cosh 9i . 

Space is taken to be periodic with period L, so each momentum is quantized: 

2rii7T 

Pi = (2-2) 

When the particle number N is large, the partition function at temperature T = 1/R and 
chemical potential fx is defined as 

< 1 > A = Z X (R, L) = tr [X N e- RHL ] , (2.3) 

where the fugacity A = e^ R . This is the standard one-dimensional free-fermion partition 
function, which is 

Z X {R 1 L)= H (l + \e- RE -), (2.4) 



where 



Thus 



. ,2tt n2 2 \ 1/2 , , 

E n = [ (— n) 2 + m 2 j . (2.5) 



F\ = —T In Z\ = -T H 1 + ^~ REn ). (2.6) 

n= — oo 

In the limit L — > oo, we can replace the sum over n with an integral, and using 

L L 

dn = — dp = — m cosh 9d6 
2tt f 2n 



yields 



mL 



J cosh61n(l + \e- e{e) )d6, (2.7) 



2ttR 

where e(d) = mR cosh 9 is the one-particle energy. 

By definition, the eigenvalues of diagonal operators vary only with particle number. 
Looking at (|2.3| ), one sees that to calculate the expectation value of such operators, one 



needs only to modify the chemical potential appropriately. Thus 



<(-!)">=<!>*=-!, (2.* 



5 



This corresponds to a chemical potential of iixT. Even though this is imaginary, ( |2.7|) is 
perfectly well-defined by analytic continuation in A. In fact, F—\ is minus the free energy 
of a free boson. The argument of (|2.7|) has a singularity at A = — 1 and 9 = 0, but the 



integral is well-defined. 

All of this analysis has been done treating the L-direction as space. When R is the 
space dimension, we recover the picture of Ql.7| ). When A = 1, Ha is the Hamiltonian 
for free fermions with antiperiodic boundary conditions. When A = (—1)^, the boundary 
conditions in the -R-direction are flipped: Ha is the Hamiltonian for free fermions with 
periodic boundary conditions. In the language of Ising spins, these correspond to periodic 
and antiperiodic boundary conditions on the spins, respectively. Thus with the equiva- 
lence evaluating < (— 1)^ > gives the ground-state energy for the Ising model with 



antiperiodic boundary conditions. 

In the limit m — > 0, this model approaches the Ising conformal field theory. To 



evaluate ( 2.7|) in this limit involves a few tricks [^3j@ which are discussed in section 3. 



Using the definition ( |TT8| ) of the ground-state energy Ea(R), the result is 

(2.9) 



El ( R ) = -i m 

u ; 2 V6.R/ 



e (-i)n(R) = [ — 



\6RJ 

The second equation should be clear even without the tricks since it is minus the energy of 
a free boson, which has c = 1. Equation ( |1.5| ) yields the scaling dimensions of the operator 
which creates these two states. Ei(R) yields c = 1/2, because h = h = 0. E(_^n(R) 
shows that the operator with conformal dimensions (1/16, 1/16) creates the lowest-energy 
state in the Ising model with antiperiodic boundary conditions. This operator is the spin 
field a. This result, of course, can be derived at the conformal point [|IU| . 

We turn to operators which interchange particles. Since there is only one particle in 
the Ising model, this is not relevant here. However, in a system consisting of two decoupled 
Ising models, there is an operator which exchanges a fermion of the first system ipi{9) with 
one from the second system ij)ii{6). This operator, which we denote by K, commutes with 
the Hamiltonian, so we can compute < K >. The Hamiltonian is diagonal on the space 
of states in theories with diagonal S-matrices, and we take the states to be orthonormal. 
Therefore, a state must be invariant under K to give a non-zero contribution to the trace 
in < K >. We define the truncated Hilbert space TCk as the space of such states; i.e., 

6 



states where for every ipj with a rapidity 9, there is also a ipu with the same rapidity 9. 
Then 

< K > x =tr HK [Af/ 2 Af/ 2 e"^" 



A A 2 , (2.10) 
where 

N N/2 

Ek = m cosh(6'i) = 2m cosh(6>i). 

i=l i=l 

The configuration sum is that of a single Ising model with N/2 particles, and the energy 
is that of a single Ising model of mass 2m. Thus 

< K > x = + \i\ne- 2E "), (2.11) 

n 

where E n is defined in ( |2.5|) . In the L — > oo limit, 

In < K > A = [ cosh6>ln(l + A/Aj/e" 2 ^)^, (2.12) 

2tvR J 



where e(9) = mR cosh(#) is still the one-particle energy. Notice that in (|2.12|) there is no 
extra factor of 2 in front of the integral. 

As before, we interchange space and Euclidean time, and interpret Ek as the ground- 
state energies for models with different boundary conditions. In the language of Ising spins, 
Hk for Xj = A// = 1 is the Hamiltonian for two Ising models, coupled only by the boundary 
condition aj(0,y) = o~u(R,y) and ai(R,y) = 07/(0, y), where ai(x,y) and au(x,y) are 
the Ising spins at site (x, y). When Xj = —1 and Xu = 1, this amounts to calculating 
< K(—1) N / 2 >. (Remember, in the truncated Hilbert space N/2 is an integer.) This 
corresponds to the boundary condition aj(0,y) = —o~u(R,y) and ai(R,y) — an(0,y)- 

In the conformal limit m — » 0, ( j2.12| ) is independent of m. Comparing this with ( |2.7| ), 
we see that 

lim (In < K > x ) = - lim F\. (2.13) 

m — >0 2 Tft — >0 

Thus we can read off the answers from ( |2 . Q| ) , and when A = 1, 

E *W - ~\ (_^) • (2.14) 
Since we have two Ising models, c = 1, and h + h = (1 — l/4)/12. Thus the field which 
creates this state has conformal dimensions (1/32, 1/32). When A = —1, 

E «w = + l (ss) ■ < 2 - 15) 

Here h + h = (1 + l/2)/12, and the conformal dimensions are (1/16, 1/16). Both of these 
results are in agreement with those found using modular transformations at the conformal 
point [f27|. 



7 



3. The TBA with a chemical potential 

In this section we briefly derive the TBA equations in the presence of a chemical 
potential, following the treatment of refs.p|fL6f. 

The starting point of the TBA is the exact two-body elastic S-matrix. For simplicity, in 
this paper we will treat only diagonal S-matrices, although the TBA can also be extended to 
models with non-diagonal S-matrices, where one must introduce massless pseudo-particles 



into the TBA system of equations [19]. The S-matrix element for the process a(6i)b(02) ~~ ► 
6(6>2)a(6>i) is denoted by S a b(0), where 6 = 61—62, the difference of the particles' rapidities. 
An S-matrix element in one spatial dimension is defined so that when two particles a{6\) 
and 6(6*2) are exchanged, the wavefunction is multiplied by S a b(6). 

For the moment, we specialize to a theory with only one kind of particle, so that there 
is only one S-matrix element S(6). The TBA applies to systems with a large number of 
particles N, where the space direction is periodic with period L. The periodicity leads to 
a quantization condition on each momentum pi = msinh(#i): 

e i ** L Y[S{8 i -0 j ) = l. (3.1) 

Defining p r {6) to be the rapdity density (i.e., p r (6)A6 is the number of particles with 
rapidity between 6 and 6 + A9), the logarithm of (|3.1|) is 



mL sinh(^) - i J d6' p r {6') In S(6 l - 6') = 27m;, (3.2) 

where an integer rii is associated with each particle. All integrals in this paper run from 
—00 to +00 unless otherwise labeled. This equation is the generalization to an interacting 
model of the one-particle relation pi = 27m;/ L. Here there is a set of coupled equations for 
the momenta. We introduce a level density p(6), so that p(6)A6 is the number of values 
of 6i which solve (|3.2| ) in the range 6 + A6. Taking the large L limit, ( |3.2| ) becomes 

p(6) = mL cosh(^) + J d6' p r (6')(j)(6 - 6'), (3.3) 

where 

.dlnS(6) 



0(61) = -i- 



86 



To determine the free energy, we use thermodynamics, as opposed to the statistical 
mechanics used in section 2. In a system with chemical potential p, (A = e^/ T ) at temper- 
ature T, we use F\ = £ — TS — pN and minimize F\ using fl3.3| ) as a constraint. Here the 
energy £ is 

£ = J Pr (6)m cosh 9d8, (3.4) 

the entropy S is 

S = J d9[phip — p r ha.p r — (p — p r ) ln(p — p r )], (3.5) 
and the particle number iV is 

N = J p r (9)d9. (3.6) 



In (|3.5| ), we have assumed that all particles are of "fermionic" type, with at most one 
particle of a given species with a given momentum 9i. This is true for all known one- 
dimensional theories, save a free boson.i Defining the pseudoenergy e{9) by@ 

Pr = PYT^> (3 ' 7) 

the minimization of F\ yields an integral equation for e(#): 

e(9) = mRcosh9 - — [ d9'(f>(9 - 9') ln(l + Ae" 6 ^). (3.8) 
2tt J 

In terms of e, the free energy is 

„ mL 



— [ d9 cosh 9ln(l + Xe-< d) ). (3.9) 

7TK J 

In the case where there is more than one species of particle, one introduces particle and 
level densities for each species a, and finds 

p a (#) = m a Lcosh(#)+^ / d9'p b r {9')<j> ah {9-9'), (3.10) 



where 



, 9 In S ab (9) 

<Pab{9) = -i . (3.11) 



3 



See [H] or [15| for a detailed explanation of this point. It is easy to treat the case where 



particles are "bosonic" — just replace p in (p.5|) with (p + p r ). 



This definition includes the chemical potential, as opposed to [16|. 
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The pseudoenergies are defined as in (|3.7|) , and minimizing the free energy yields the 
integral equation 



e o (0) = m a R coshfl - ^ i- / d0'<M0 - 0') ln(l + A 6 e 



-6 b (fl)^ 



(3-12) 



We define the scaled free energy density f\(mR) as 



f x (mR) = ^F A , 



and it is given by 



/ A (m J R) = ^/ A °(m J R), 



f%(mR) = - — [ d6 cosh 6 \n(l + \ a e- £ ^ e) * 
2n J 



(3.13) 



The conformal limit of (|3.13|) is the limit mR — > 0. To evaluate this, we define 
e fcmfc (#) ee e(9 — x), where x = ln(m-R/2). When x — > — oo, this yields equations which do 
not depend on mR: 



ef nk {6) =e e -J2^J d6'<p ab (6 - 9') ln(l + \ h e 



'(0) 



and 



7T 



/a(0) = — / e 9 ln(l + A a e- 



)■ 



(3.14) 



(3.15) 



There are a few tricks which enable us to evaluate this without knowing the full function 
e(9). We take the derivative of (|3.14|) and substitute for e d in (|3.15|) . The (f) a b can be 
removed from the resulting expression by using ( |3.14[ ) again. Every term depends on 9 only 
through e ktnk , and also has a de ktnk /d9 in it. Since the shifted e(9) is always decreasing, 
the integral over 9 can be replaced by an integral over e kink , yielding 



2/a(0) 



r, u W 



eX a e~ 



.1 + A a e 

where the contour C a runs in the complex-e plane from e k 



+ ln(l + \ a e 



-oo) to e 



kink . 



(3.16) 



oo). In theories 



without massless particles, e kink (oo) = oo. To determine e /cinfc (— oo), we use the property 



kink i 



of (|3.14| ) that e fcinfc (6») is flat from 9 = -oo to around 9 = 0, and that (p a b falls off 
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exponentially. This enables us to pull ln(l + A&e 6b (°°^) out of the integral in ( |3.12|) , 
which yields a system of equations for x a = exp(e^ mfc (— oo)) = exp(e a (0)): 



x b 



where 



N ah = - 



1 

2r7 



(3.17) 



(3.18) 



When A a = 1, the contour C a lies on the real axis, but in the cases with imaginary chemical 
potential we will discuss later, it may not. Of course, one can deform C a , paying attention 
to the logarithmic cuts in ( |3.16| ) starting at /i a — e a = i(2n + l)ir. 
We write ( |3.16|) in the form 



A(o) 



TV * — ' 



(3.19) 



When A a are either 1 or —1, and all x a > 1 for a where A a = — 1 , then the contours C a 
run along the real axis and do not go through the logarithmic singularity. In this case, a 
change of variables in the integral ( ^.16| ) yields 



A 0*0 =£(1/(1 + x)) 
C-x(x)=-L(l/x), 



(3.20) 



where 



L(x) 



dy 



Iny ln(l-y) 



+ 



y 

There are many known identities involving the Rogers dilogarithm function L(x), and they 
can often be used to evaluate ( |3.19| ) exactly ||T3]| p8|. 

Since this derivation allows an arbitrary chemical potential for each particle, the ex- 
pectation value of any symmetry operator diagonal on the space of particles is obtained 
by setting A a to be the value of A on the particle a. A variety of such examples will be 
presented in sections 5 and 6. 



4. Non-diagonal operators 

In this section, we show how to calculate the expectation value of a symmetry operator 
K which is not diagonal on the space of particles. This involves truncating the Hilbert 
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space, just like in the (Ising) 2 model in the second half of section 2. The procedure for 
calculating 

< K >= tr [Ke~ RHL ] (4.1) 

basically follows that of section 3: we find the appropriate quantization of momenta in 
periodic boundary conditions, and use this as a constraint while minimizing the free energy. 
This results in an integral equation for the pseudoenergy e(6) and an expression for the 
free energy in terms of e(6). If there are no particles neutral under K, the asymptotics 
of these equations are obtainable, allowing the calculation of the conformal dimensions of 
the operator which creates the excited state. In others, it seems that one must resort to a 
numerical calculation. 

The relation which determines the level densities p a (9) remains (3TTU) , but we must 



modify the free energy. The Hamiltonian is diagonal in theories with diagonal S-matrices, 
so only states invariant under K contribute to the sum in Thus, as before, < K > 

is the sum of e~ RHh over states on which K = 1. This restriction changes the entropy. 
The species of particles form multiplets under K, and all species of particles in a multiplet 
must have the same level density and particle density (this is true even in the unrestricted 
case). The crucial effect of the restriction is that the entropy from only one species of 
each multiplet contributes to the free energy. This is because in order to specify a given 
state with K = 1, one needs only to specify the set of momenta for one species in a given 
multiplet; the other sets must be identical by K invariance. Minimizing the free energy, 
one finds the coupled integral equations are modified to 

e a (6) = n a m a R cosh 9 -J^^ f d9'cj> ab {e - 6') ln(l + A^e"^), (4.2) 



where n a is the number of species in the multiplet a. The scaled free energy, defined as 
R 

L 



f K {mR) = ^ln<K > x = RE K (R), is 



(4.3) 

_ _ mR r ddcQshdln ^ + Aa e- e aW). 

27T J 



We do not treat operators which take a particle into a superposition of particle states. 
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In the limit mR — > 0, doing the same substitutions as in in section 3, one finds that 
(specializing to the case where all n a are 1 or 2) 

/*(0) = -c A (£)-!l>' (4.4) 

where the sum over A is only over particles with ua = 2, and c\ is the result obtained 
without truncating the Hilbert space; e.g., c\ is the central charge derived by the methods 
of section 3 with A = 1. Unfortunately, I can find a closed- form expression for /ir(0) like 
( |3.16| ) only in the case where all the n a are the same. In the case where all n a = 2, (|4.4| ) 
yields 

/*«» = («) 

When A = 1, (P-.5Q gives h + h = (|c)/12, and the operator which creates this excited 
state has conformal dimensions (c/32,c/32). This is quite an odd result — it applies to 
any theory with a conjugation symmetry that leaves no particles neutral. This includes 
the situation where two copies of a model are coupled by boundary conditions, as in the 
(Ising) 2 model treated in section 2. 



5. Examples from Toda minimal S-matrices 

5.1. The Three-State Potts Model 

The three-state Potts model at its conformal point has central charge c = 4/5. Per- 
turbation by the thermal operator leaves the theory integrable, and its exact S-matrix 
was conjectured in [p9fl . The particle spectrum consists of a particle (labeled 1) and its 
antiparticle 1. These particles are in the two-dimensional representation of the model's 
S3 symmetry. This symmetry is generated by the operators A, which is diagonal and 
multiplies particle 1 (1) by e 27 ™/ 3 (e -27 ™/ 3 ), and K, which exchanges the two particles. 

Calculating < K > requires merely reading off the answer from section 4. The con- 
formal dimensions of the operator creating the associated excited state are found from 
equation (|4.5| ), and in this case are (1/40,1/40). In the language of the lattice Potts 
model (where a spin a takes the values 1, e 27 ™/ 3 or e -2771 / 3 ), this corresponds to find- 
ing the lowest eigenvalue of the Hamiltonian with "twisted" boundary conditions, where 
cr(0, y) = o~*(R, y). This is in agreement with the conformal result |10]] ■ We note that even 
though the (1/40,1/40) operator is not part of the operator content of the three-state Potts 
model on the torus, it is present with these boundary conditions. This is analogous to the 
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situation in the Ising model, where the free fermion is not part of the toroidal operator 
content, but does belong in the case of antiperiodic boundary conditions |10 . 



Calculating < A > requires the introduction of a chemical potential, as in section 3. 
In this case, we must have a complex fugacity: Ai = e 27 ™/ 3 and Ay = e -27 ™/ 3 . We define 
the TBA equations here by analytic continuation from Ai = Ay = 1, whose TBA system 
was derived in ||. In this analytic continuation, we keep Ai = (Ay)*, enabling us to require 
that e 1 = (e 1 )*. This specifies which branch of the logarithm is taken in ( |3. 12|) and ( |3.13|) . 
The free energy ( |3.13| ) remains real because fl = (fx*)*- These resulting TBA equations 
are those discussed in ||. 

To evaluate / in the conformal limit, we must find x a by solving ( [3.171) . In this case, 
Nn = N22 = 1/3, and N12 = N21 = 2/3. || This equation has more than one solution, but 
the correct one is the continuation of the A = 1 result of xi = xj = (y/b + l)/2. It is 



Xi Xj 1 



-7T2 „— §7Tj 2 



which with our analytic continuation yields 



.(1 -V5), 



I (5-1) 
x T =-(VE- l)e + *t. 



Thus the contour C\ in (|3.16| ) starts at \nx\ = ln[(-\/5 — l)/2] — in/3 and goes off to 



positive infinity, asymptotically approaching the real e-axis from below. It goes over the 
logarithmic singularity at e = in/3. The contour Cy is of course the complex conjugate of 
C\. The actual shape of these contours is determined by the solutions e a (9), but we do 
not need this information to evaluate ( |3.16| ). We can deform the contour C\ so that Im(e) 
is always —in/3. This goes right through the singularity of the integrand, but the integral 
is well-defined. Shifting e — > e + in/3 in ( |3.16|) yields 



2/a(0) = -- T de f-ill + ln(l - e" e ) + 

71 ./in |a?i I L 1 - e V 3 / 



In \ xi I 

+ complex conjugate. 



(5.2) 



The last term in the square brackets is a result of the shift, and can be evaluated explicitly, 
giving 

I 1 \ 7.7T 

+ complex conjugate. 



1 \ in 

- T lnd-e- 



In \x\ I 
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This is not zero because the argument of the logarithm at e = In |xi| is negative, and the 
way we have analytically continued means that the logarithm from C\ gives — in, whereas 
the logarithm from Cj gives +irr. Thus the total contribution of this term is — 27r/3. To 
evaluate the other two terms, we split the contour into two parts, one running from In |xi| 
to zero, and the other from zero to infinity. Then the change of variables x = 1 — e e in the 
first part and x = e~ € in the second yields 



Mo) = -- 

7T 



TT 2 



-2L(i(3-v / 5))-2L(l)+ ; , 



(5.3) 



Since L(l) = tt 2 /6 and L(±(3 - V$)) = tt 2 /15, this yields 



E * - + l (sg) • < 5 ' 4 > 

The operator which creates this state has dimensions (1/15, 1/15), as was derived in ||. In 
the lattice Potts model, this corresponds to "cyclic" boundary conditions, where cr(0, y) = 



exp(27u/3)or(.R, y). This result has been derived at the conformal point [10 



5.2. Z n _|_i parafermions 

The Ising model and the 3-state Potts model are the first two models in the hierarchy 
of Z n+ i parafermions |[30|| . At the conformal point these models have central charge c = 



2n/(n + 3). There are "spin" fields (k = l...n) which have conformal dimension 
h = h = fc(n + l-fc)/2(n + l)(n + 3),and "C-disorder" fields $ s (s = 0, 1, 2 . . . < (n + l)/2), 
which have conformal dimensions h = h = [n — l + (n+l — 2s) 2 ]/16(n + 3). These conformal 
field theories remain integrable when perturbed by the primary field with dimensions 
(2/(n + 3),2/(n + 3)), and the S-matrix is conjectured to be the A n minimal Toda S- 
matrix|3l]].i There are n particles with mass m a = sin(7ra/(n + 1)), a = 1 . . . n. 

The value of iV a b for a simply-laced Toda minimal S-matrix are given by the matrix 
relation iV = 1(2 — I), where / is the incidence matrix (= 2 — C, where C is the Cartan 
matrix) of the Lie algebra corresponding to that Toda theory. It was shown in |33|] that 



the relation ( |3.17| ) for simply-laced Toda theories can be transformed to the simpler form 

X a = H(\a + X a ) 1 z I « b . (5.5) 



For arguments as to why this is so, see [32]. 
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The symmetry of a Toda S-matrix is the symmetry of the extended Dynkin diagram. 
For A n this is the dihedral group P n +i 5 and is generated by two operators A and K. The 
diagonal Z n+ i symmetry is generated by A, which when acting on the particle a gives the 
phase exp(2-7ua/ (n + 1)). K implements a Z2 symmetry, which exchanges particle a with 
the particle n + 1 — a. For more details on these and other Toda S-matrices, and on the 
TBA in these theories, see [ITH] and |]33||. 



By using the method of section 3, we can determine < A k > by setting the chemical 
potentials A a = exp(z27r/ca/(n + 1)). For the A n Toda minimal S-matrices, the solution of 

(H) is 

x S in ih±^IL sin M»I 

•^a _ n+3 n+3 /r n\ 

A a Sin pj— 

n+3 

These values should be used with (|3.16| ) to determine the ground-state energy in the 
conformal limit. In the Ising (n=l) and three-state Potts (n=2) models, the excited state 
associted with < A h > is created by cr fc . Therefore, we conjecture that this holds true for 
this entire hierarchy. This has been checked for n = 4 in || . In the next section, the n = 3 
case is discussed in detail. However, I was unable to find a general way of evaluating the 
integral in ( |3.16[ ), and evaluating them case- by-case is straightforward but tedious. 



By using the method of section 4, we find the integral equations which determine 
< K >. For n even, there are no particles neutral under K, so we can read off the answer 
in the conformal limit from (|4.5| ). The operator which creates the associated excited state 
has dimensions n/16(n+3), and is the $ n /2 field. To evaluate < K > in the conformal limit 
for n odd would require solving the integral equation ( |4.2| ) numerically, and substituting 
the result into (fO|). However, it can be easily estimated using 

M0)^4 L (TTM (5-7) 

where L is defined in ( |3.20| ), and the x a are those of the untruncated model. (The x a do 
not change because ( |3.12| ) and ( f4~2^ ) are identical in the m — > limit.) Evaluating the 
dilogarithms numerically, we find that 

MO) » -\ (I 

for all odd n. The associated field has dimensions (n — l)/16(n + 3), and is $( n +i)/2- One 
should be able to verify this in the conformal limit by using the analysis of ||34|| . 
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5.3. 7i± parafermions (in detail) 

In this section we calculate the conformal limit of < A > and < A 2 > for Z4 
parafermions. The S-matrix is that of the A3 Toda theory, where the mass spectrum 
consists of three particles 1, 2 and 3, with masses 1/a/2, 1 and respectively. 

To calculate < A 2 >, we set A2 = 1, Ai = A3 = —1. The solution of (|5.5|) is given in 
(|5.6|) , and is x\ = X3 = 1, and x^ = 0. It is not the unique solution, but it is easy to verify 
that this in fact is the solution obtained by deforming Ai and A3 continuously from 1 to 
— 1. Since all the A a are either 1 or —1, we use (|3.19f) and ( |3.20| ) to obtain 

/a*(0) = --[L(1)-2L(1)] 

(5-8) 



.6. 

since L(l) = tt 2 /6. The equivalence flL|) means that this can be interpreted as the lowest 
eigenvalue of the Hamiltonian with boundary conditions appropriate to A. (The lattice- 
model interpretation of these boundary conditions will be given in the next subsection.) 
The conformal dimensions of the operator creating this state are given by the asymptotic 
formula ( |1.5|) , and are (1/12,1/12). This operator is the spin field 02. 

The calculation of < A > in the conformal limit is a bit more involved. We set 
A a = (i) a , and the solution of ( |5.5p is X\ = 0, x% = 1. Using ( |3.20| ), we find 

7T 6 

Since x± and X3 are real the contours in ( |3.16| ) run along the real axis: 



C Xl +£\ 3 = -— I dc 

Ztv 



3C 



+ ln(l + ie 



+ complex conjugate. (5.9) 



1 + ie 

Even though the two terms in the square brackets individually diverge when e — ► 00, their 
sum does not, and the integral is well-defined. Since the integral is over the real axis, we 
can add it and its complex conjugate, obtaining 



1 r°° r 2e 



C Xl +C X3 = -— I de 



00 



ee , ^ _ 2ex 



l + e- 



-2e 



+ ln(l + e- 2e ) 



Ztv 



(5.10) 



where we do the change of variables x = 1/(1 + exp(2e)) to obtain the dilogarithm in the 
second line. Adding all the contributions together, we obtain 

fA(0) = \ (I) , (5.11) 

and associated operator is the spin field crj, with conformal dimensions of (1/16,1/16). 
We also have < A 3 >=< A >, which agrees with the general parafermion symmetry 
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5-4- D n Toda theories 

The D n Toda S-matrices describe a hierarchy of models which in the conformal limit 
all have c = 1. They correspond to a scalar field with radius \Jn/2 orbifolded by its Z2 
symmetry []35[]. They are the continuum limit of the lattice Ashkin- Teller model, which 



consists of two Ising spins at each site with a four-spin coupling specified by the value of 
n. The (Ising) 2 model and the Z4 parafermions are the first two models in this series, 
since D2 = A\ + Ay and -D3 = A3. The D n model has n particles, with mo = = 1 and 
m a = 2 sm(na/2(n — 1)), where a runs from 1 to (n — 2). As before, the symmetry is that 



of the extended Dynkin diagram. For more details of these models, see |36| or |IS]. The 
TBA was applied to these models in [fL5| . 

The symmetry structure is slightly different for n odd or even. For all n, there is 
a Z 2 generated by K, which interchanges with 0. (For n=4, this becomes £3, which 
interchanges 1,0 and 0.) For n even, the diagonal symmetry is Z2 x Z2, and is generated 
by two operators (A±, A2), which act on the particles with eigenvalues a\ and «2, where 
a 1 (a) = 02(a) = (— l) a , cki(0) = 02(0) = 1, and «i(0) = 02(0) = —1. For n odd, the 
diagonal symmetry is Z 4 , and is generated by A, with eigenvalues a (a) = (— l) a , a(0) = i 
and a(0) = —i. 

To calculate < A1A2 > for n even or < (A) 2 > for n odd, we set Aq = Aq = —1 and 



A a = 1. The solution of (|3.17| ) is 



X = Xq = 1 

sin ((£ri)ZL )s i n( ^±i)2E) (5.12) 



sin 2 (^) 



Since all A a are ±1, ( |3.19|) and ( p.20| ) are applicable in the conformal limit. Notice that 
x'i = and that the x a for a > 2 are those for the A n _ 3 Toda theory, as displayed in (|5"T6| ) 
with k = 0. Thus 

w»)=-;i M +i-i-i) 

6 " (5.13) 
= -|(l--)- 

The associated excited state is created by the field with conformal dimensions (l/4n, I/An). 
This corresponds to anti-periodic boundary conditions on both Ashkin- Teller spins, and 
has been derived in the conformal limit [BTj . 
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To calculate < A\ >=< A 2 > for n even, we set Ao = 1, Aq = —1 and A a = ( — l) a . 
The solution of (|3T7f ) is 

x a =1 a odd (5.14) 



x a =0 a even. 



Thus (FugD and (ET20D yield 



/ Al (0)=-^(-l + i + l-l + ... + l-l) 

, 5 2 (5.15) 

1 /7T X 



2 V6 

The associated excited state is created by the field with dimensions (1/16,1/16) and is 
called the spin field. This corresponds to anti-periodic boundary conditions on one of the 



Ashkin- Teller spins, and has also been derived in the conformal limit [27]. 

Calculating < A > for n odd requires setting Ao = i, Aq = — i and A a = (— l) a . The 
solution of ( |3.17| ) is ( |5.14j ) in this case as well, but we cannot use ( |3.19| ) this time. However, 
evaluating (|3.16|) for the three particles 0, and 1 is identical to the calculation of < A > 
done for the D3 model (Z4 parafermions) treated in the previous subsection. The rest of 
the terms in (|3.16| ) can be evaluated with ( |3.19| ), so we find 



A(0) = -£(-i + i-i + ... + i-i) 

, 6 2 (5.16) 

1 /7T X 



2 V6 

and the associated field has conformal dimensions (1/16, 1/16) here as well. 

Evaluating < K > exactly for n > 2 would require a numerical calculation. However, 
using the approximate formula ( |5.7|) yields conformal dimensions of 



fh h)&( —L(l/n), — L(l/n) ] w (— , — ), (5.17) 
which is the result obtained at the conformal point. 

5.5. Eg, Ej Toda theories 

The tricritical Ising model (c = 7/10) perturbed by its thermal operator is believed to 
be described by the seven-particle £7 Toda S-matrix|]37]]. The solution of (|5.5| ) for A a = 1 
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is given in [15|. This model has a diagonal Z2 symmetry which multiplies particles 1, 3 
and 6 (in the notation of |T5| ) by —1. The solution of ( |5.5|) for Ai = A3 = A6 = — 1 is 



£4 = x-j =0 

xi = x 3 = x & =1 (5.18) 
1 + ^5 

X2 = X 5 = . 



Thus CT ) and (Q) yield 



/a(0) = -^-f-l + |- l + l + |- l + l 

Z7T \ 5 5 

1 /7T 

~5 V6 



(5.19) 



which means that the associated operator has conformal dimensions (3/80,3/80). 

The tricritical three-state Potts model (c = 6/7) perturbed by its thermal operator is 
believed to be described by the six-particle Eq Toda S-matrix|38[]. The solution of ( |5~5| ) 
for A a = 1 is given in This model has an S3 symmetry generated by the operators A 
and K. The diagonal operator A multiplies particles 1 and 3 by exp(27ri/3), 1 and 3 by 
exp(— 2ni/3), and leaves particles 2 and 4 neutral (using the labels of Hl5"j). K interchanges 
1 with 1 and 3 with 3, while leaving 2 and 4 neutral. The calculation of /a(0) was outlined 
in ||, and the associated operator has dimensions (1/21, 1/21). To calculate /k(0) would 
require a numerical calculation, since there are particles left neutral by this symmetry. 
The approximate formula (|5.7| ) yields the conformal dimensions of (1/56,1/56) for the 
associated operator. 



6. Examples from the minimal models 

6.1. The massless perturbation of the tricritical Ising model 

The conformal tricritical Ising model (the n = 4 minimal model) perturbed by the 
0i,3 operator 

Stci — > Stci + 9 J d 2 x(pi j3 

is described by two different field theories, depending on the sign of g. For g positive, one 
obtains a theory of massive kinks with a non-diagonal S-matrix||39||. Since the S-matrix is 



not diagonal, one must introduce pseudoparticles into the TBA[19]. Due to the difficulty in 
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undersanding the effect of boundary conditions on the pseudoparticles, we will not discuss 
this case here. When g is positive, the model flows to the Ising model (n = 3 minimal 
model). Thus there must be massless excitations, which are the only particles which will 
remain in the infrared limit. As shown in |20|1 , one would expect there to be a right- moving 
and a left-moving massless fermion, which in the infrared limit become the free Majorana 
fermion of the Ising model. The S-matrix for these fermions was conjectured in ||2U[| , and 
the TBA system was derived. Since these particles are massless, the rapidity relations 
( pTT|) do not apply but instead are replaced by 

M e 

f right 

M , ™ 
Pleft = -y« . 

where M is the scale of the theory, which depends on g. The TBA equations ( |3.12| ) are 
changed accordingly. 

We calculated < (—1)^ > in the Ising model in section 2, and showed that it corre- 
sponded to periodic boundary conditions in the indirection on the fermion. The same is 
true here. The calculation can be done by introducing fugacities X r i g ht = ^left = — 1- The 
analysis of section 3 is easily repeated to take into account the modified rapidities(|Q|). One 
then obtains the TBA equations conjectured in || and as describing the excited-state 
energy for this model. We thus have provided a physical interpretation of these equa- 
tions. These papers show that the operator which creates the excited state has dimensions 
(3/80,3/80) in the conformal limit. 

6.2. The "staircase" model 

A simple TBA system with many intriguing similarities to the minimal models has 



recently been studied |21|]. The system consists of a single particle $ with mass m and with 
S-matrix 

s(e) = sin ll:' C03h h e ; , (6.2) 

smh t) + i cosh t>o 



where 9q is a free parameter. Using (|3.11|) and (|3.18|) , one finds = — 1. Solving ( |3.17p 



yields x$ = 0, which means that in the m — > limit, c = 1. This S-matrix is the analytic 



continuation of the S-matrix for the sinh-Gordon model ||4C 

1 



shG 



(<9 M $) 2 - 2^cosh/3$ 
21 



d 2 x. (6.3) 



to complex coupling 7 = £ + iOo, where 

/3 2 /8 
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1 + /3 2 /8tt' 



In []21| it is shown that at large 6> , f(mR) develops an unusual "staircase" structure. 
For all other known models, the structure of f(mR) (as a function of mR on a log scale) 
is a plateau at /(0), which eventually falls off smoothly to f(oo) = 0. One then can read 
off the central charge of the conformal theory from the value of f(0). In this "staircase" 
model at large 6>o, f(mR) develops plateaus in various regimes of mR, making the plot 
of f{mR) look something like a staircase. Even more intriguing is that the values of 
f(mR) on the plateaus are the values which give the central charge of the minimal models, 
c = 1 — 6/n(n + 1). Thus the picture one has is that one starts out at a c = 1 conformal 
theory at m = 0, and then as m is increased, the renormalization group flow takes one 
into the neighborhood of the fixed point corresponding to each of the minimal models. 
The trajectory spends some time in the neighborhood of each model, before the instability 
takes over, forcing it to flow to the next model. This goes on until it finally reaches the 
c = trivial model, where all the degrees of freedom have been scaled away. The staircase 
structure was derived numerically and analytically. It is derived analytically by showing 
that the TBA system defined by ( |6.2| ) can be accurately approximated by the TBA system 
conjectured in |T9[ as describing the flow of the ground state from the n-th to the (n — l)-th 
minimal model[^TJ. For a much more detailed description of the staircase model, see the 
original paper pTj] . 

Since there is no pole in the S-matrix ( |6.2j ) corresponding to $ as a bound state of 
itself, this model must have a Z2 symmetry <f> — > — $. This is also apparent from the 
sinh-Gordon description ( |6.3| ). Thus in our usual manner we can set A = — 1 and calculate 
< (-1) N >. The solution of ( ftTFD is = 2, so 

/-i(0) = ^)- (6.4) 

Thus the associated operator in the conformal limit has dimensions (1/16, 1/16). It is easy 
to repeat the argument of [^T| to show that fi(mR) develops plateaus, only this time the 
TBA system can be approximated by the TBA system studied in sections 5.1 and 5.3 of 
ref. 0. This is identical to the system conjectured for the ground-state energies JH5J, except 
in this case all the A a = — 1 instead of 1. The plateaus do not form a staircase here; 
for n odd (even) the plateaus are at values greater (less) than /_i(0). The plateaus take 
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values of f-i corresponding to the conformal dimensions of the field <fi n /2,n/2 f° r n even 
and 0( n +i)/2,(n+i)/2 f° r n odd. This is conjectured to describe the flow from the n-th to 
the (n — l)-th minimal model of an excited-state energy for n even, but not n odd||. This 
is for several reasons. The first is that the numerical results for n odd beyond lowest order 
do not agree with perturbation theory around the conformal point. The second is that one 
expects the flows in minimal models to take the <f>n operator in the n-th minimal model 



to the operator in the [n — l)-th minimal model|4^] |41]; this is true here only for n 
even. 

So while our result shows that excited states in the staircase model resemble those in 
the minimal models in some respects, there is a crucial difference. Although the flow down 
the chain of operators seems very natural from the point-of-view of the staircase model, 
its relation to flows in the minimal models is still not completely clear. 



7. Conclusions and speculations 

We have seen that the TBA is useful for calculating not only the Casimir energy of 
the lowest-energy state in the theory, but for a number of excited states as well. This is 
done by choosing the boundary conditions to project out the ground state. In the limit 
R — > oo, all of the states discussed here have E — > 0, and become degenerate with the 
ground state.@ It is not clear if it will be possible to project out all such states and obtain 
information on massive states, in the manner of 03] [13] . 

The most obvious open direction is to understand how these results can be used to 
study excited states in theories with non-diagonal S-matrices, where one must introduce 



massless pseudo-particles into the TBAjjT^]. TBA equations in some such situations have 
been conjectured, but it is not obvious how to interpret these equations as boundary 
conditions on the field theory. However, this hopefully is only an interpretational issue: 
all the formalism should be applicable to these cases. 

There are many other operator expectation values that one can calculate. For example, 
it is straightforward to calculate < (—l) Na > for any particle a in a diagonal scattering 
theory, since (—l) Na is a symmetry of the S-matrix even if it is not a symmetry of the field 
theory or of the fusion rules. The interpretation of this in terms of boundary conditions 
is not immediately obvious. However, in the minimal models on the lattice (on and off 



7 I thank Tim Klassen for pointing this out. 
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criticality) it is possible to define boundary conditions so that any of the primary operators 
creates the lowest-energy state[^5j. This seems to be true for all rational conformal field 



theories as well [|TT| . It is thus not difficult to believe that once suitable ways are found 
of implementing the boundary conditions, one will be able to calculate the energy for all 
such excited states by calculating operator expectation values with the TBA. 
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